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Abstract — In this paper, we propose a selection region based 
multihop routing protocol with directional antennas for wireless 
mobile ad hoc networks, where the selection region is defined by 
two parameters: a reference distance and the beam width of the 
directional antenna. At each hop, we choose the nearest node to 
the transmitter within the selection region as the next hop relay. 
By maximizing the expected density of progress, we present an 
upper bound for the optimum reference distance and derive the 
relationship between the optimum reference distance and the 
optimum transmission probability. Compared with the results 
with routing strategy using omnidirectional antennas in fl), we 
find interestingly that the optimum transmission probability is a 
constant independent of the beamwidth, the expected density of 
progress with the new routing strategy is increased significantly, 
and the computational complexity involved in the relay selection 
is also greatly reduced. El 

I. Introduction 

Determining the capacity region of wireless networks has 
been an open problem for more than half a decade. In the 
seminal work Q, Gupta and Kumar proved that the transport 
capacity for wireless ad hoc networks, defined as the bit-meters 
pumped every second over a unit area of the network, scales 
as G(y/n) in an arbitrary network, where n is node density. 
In 0, Weber et al derived the upper and lower bounds on 
the transmission capacity of spread- spectrum wireless ad hoc 
networks, where the transmission capacity is defined as the 
product between the maximum density of successful trans- 
missions and the corresponding data rate, under a constraint 
on the outage probability. All these works address the single 
hop transmission. Recently, Baccelli et al |4] proposed a 
spatial reuse based multihop routing protocol, and derived the 
optimum transmission probability. In their protocol, at each 
hop the transmitter selects the best relay to maximize the 
spatial density of progress. By assuming each transmitter has 
a sufficient backlog of packets, Weber et al in [5 ] proposed 
the longest-edge based routing protocol where each transmitter 
selects a relay that makes the transmission edge longest. In (6), 
Andrews et al defined the random access transport capacity 
with the maximum allowable total number of transmissions 
per packet along multiple hops. In d, Li et al proposed a 
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selection region based multihop routing protocol to guarantee 
the message transmitted towards the final destination, where 
the selection region is defined by two parameters: a selection 
angle and a reference distance. By maximizing the expected 
density of progress, the author derived an upper bound on the 
optimum reference distance, and the relationship between the 
optimum reference distance and the optimum selection angle. 

The above literatures only considered the wireless networks 
with omnidirectional antennas. In Q, Yi et al investigated 
the wireless networks' capacity using directional antennas, 
extending the network capacity with omnidirectional antennas 
in (2l to that with directional antennas. In (8), Spyropoulos 
et al discussed the network capacity gain one can achieve by 
using directional antennas over that by using omnidirectional 
antennas and how these bounds are affected by important 
antenna parameters like gain and beamwidth. In [9], Dai et 
al combined the multiple channels and directional antennas 
together and shown that they improve the network capacity 
due to the increased network connectivity and reduced inter- 
ference. However, previous literature usually focuses on the 
scaling laws of the network capacity. 

In this paper, we extend the former work in [ 1 ] with omnidi- 
rectional antennas to that with directional antennas, and derive 
the close-formed expected density of progress of the network, 
which is defined as the number of packets progress toward 
their destinations in a unit area of the network. Compared with 
the routing strategy using omnidirectional antennas in [1], due 
to the directional antennas, the selection region based routing 
is implemented much easier and the calculation burden for 
nodes to select the relay is also decreased. 

The rest of the paper is organized as follows. Directional 
antenna model, network model and routing strategy are de- 
scribed in Section II. The optimization for selection region 
and transmission probability are presented in Section III. 
Numerical results and interpretations are given in Section IV. 
Finally, Section V summarizes our conclusions. 

II. Network Model and Routing Protocol 

In this section, we first present the simplified directional an- 
tenna model, then define the network model and the selection 
region based routing protocol using directional antennas. 



A. Directional Antenna Model 

In the study of wireless networks, the antenna model is often 
grouped into omnidirectional and directional. Omnidirectional 
antenna radiates signals equally well in all directions, while 
directional antenna has gain in the direction of the main lobe 
at which it is pointing. Thus, with directional antennas the 
interference can be decreased and nodes located in each others 
neighborhood may transmit simultaneously, which increase 
spatial reuse of the channel. To simplify the analysis, we model 
the power pattern of the directional antenna as a circular sector 
with angle (p, where (p is the beamwidth of the antenna, see 
Fig. 1. In the following analysis, we assume the transmitters 
with directional antenna and the receivers with omnidirectional 
antenna, which is called Directional Transmission and Omni- 
directional Reception (DTOR) as mentioned in (7). 

B. Network Model 

Assume nodes in the network follow a homogenous Poisson 
Point Process (PPP) with density A, and slotted ALOHA as 
the medium access control (MAC) protocol. During each time 
slot a node chooses to transmit data with probability p, and 
to receive with probability 1 - p. Therefore, at a certain time 
instant, the transmitters follow a homogeneous PPP (II t x ) with 
density pA, while the receivers follow another homogenous 
PPP (T1r x ) with density (1 - p)A. At each hop in the multihop 
transmissions, a transmitter tries to find a receiver in U Rx as 
relay. We consider the nodes are mobile, to eliminate the 
spatial correlation, which is also discussed in fl4). We also 
assume that all transmitters use a fixed transmission power 
p and the wireless channel combines the large-scale path- 
loss and small-scale Rayleigh fading. The normalized channel 
power gain is given by 



G(d) = ^-, 
d a 



(1) 



where y denotes the small-scale fading, drawn from an expo- 
nential distribution of mean ^ with probability density function 
(PDF) f y (x) - p exp(-yux), and a > 2 is the path-loss exponent. 

For the transmission from transmitter i to receiver j, an 
attempted transmission is successful if the received signal-to- 
interference-plus-noise ratio (SINR) at the receiver j is above a 
threshold J3. Thus the successful transmission probability over 
this hop with distance d\j is given by 



Pr 



(2) 



Xkeu TX \{i}P7idkj a + 11 

where i e Il Tx , j e U Rx . Since we use directional antenna 
with beamwidth (p for transmission, the interfers seen by 
a specific receiver follow a homogeneous PPP (Ht X ') with 
density p§^A. Thus Tjken Tx ,\{i} P7idkf a is the sum interference 
seen at the receiver j, where dkj is the distance from interferer 
k to receiver y, and n is the average power of ambient 
thermal noise. In the sequel we approximate n = 0, which 
is reasonable in interference-limited ad hoc networks. From 
L4j, the successful transmission probability from transmitter i 
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Fig. 1. Selection region with the directional antenna 



to receiver j is derived as 

P s = exp[p^Atdij 2 y 

27T 2 /a 



where 



t = 



2/a 



sin(27r/of) 

C. Routing Strategy with Directional Antennas 



(3) 
(3a) 



Considering a typical multihop transmission scenario, where 
a data source sends information to its final destination that 
is located far away, and it is impossible to complete this 
operation over a single hop. Thus a multihop transmission is 
needed. In multihop wireless networks, if we assume position- 
determined relays exist to ensure each hop shares the same 
distance that aggregates to form the path from the data source 
to its final destination, the optimum transmission distance at 
each hop is derived in (6). In this case, the transmission 
distance is used to determine the location of a relay. In a 
practical case, nodes are usually randomly distributed, thus 
relays may not be located just over the optimum transmission 
distance. To guarantee a relay existing at a proper position, we 
use the selection region based multihop routing protocol with 
directional antennas. For each transmitter along the route to the 
final destination, the selection region includes two parameters: 
the beamwidth (p and a reference distance r m , as shown in 
Fig. 1, where the selection region is defined as the region 
that is located within angle (p and outside of the distance 
r m . Here, the transmitter is located in the circle center O, 
IBOC = lAOC = (p/2, and OC points to the direction of 
the final destination. At each hop, the transmitter selects the 
nearest receiver located in the selection region as the relay. 

Compared with the routing strategy using omnidirectional 
antennas in [1|, the selection region based routing with direc- 
tional antennas can be implemented much easier. Since the 
transmitters are equipped with directional antennas, only the 
receivers within the angle (p can receive the radiated signals 
from the transmitters, thus the process to delimit the potential 
receivers within the angle (p do not need any calculation. 
However, for the routing protocol with omnidirectional an- 
tennas, to determine if a receiver is located within the angle 
(p or not needs a complicated calculation, e.g., calculating the 
angle (p between the line from the transmitter towards the final 
destination and the edge from the transmitter and the potential 
receivers, and making comparisons between these cp with (p to 



decide which nodes are located within the selection region. 



III. Reference Distance and Transmission Probability 
Optimization 

In this section, we derive the optimum values of the 
transmission probability p and the reference distance r m for 
different beamwidth ip by maximizing the expected density of 
progress. 



A. Upper Bound for Optimum Reference Distance 

As in (4), the density of progress is defined as 

D = pA • P s • dcos0, 

where P s is the successful transmission probability defined 
in (2), d cos cp is the projection of the transmission distance 
d along the line connecting the transmitter and the final 
destination. Since the receivers follow a homogeneous PPP 
with density A(l - p), the cumulative distribution function of 
the transmission distance d is given as 



Pr(d < r) = 1 - exp 



-A(\-p)^{r 2 -rl) 



,r m < r < oo. (4) 



Since is uniformly distributed over [—<p/2,<p/2], which is 
independent of d, the expected density of progress is given by 

<p 

CO 2 

E[D]=pA 2 J J e- p % Atx2 x cos (f>f d (x)d(f>dx 

r _ <P 

m 2 

= A 2 p(l -p)r||^Jr 3/2 exp(i(l -p)| r 2)sin(|), 

(5) 

where fd{x) is the probability density function of d obtained 
from (4), k = % (f + (1 - p)), t is defined in (3a), and 

oo 

r(|,£r^j = J e~ f xi~ l dx is the incomplete Gamma function. 

kr 2 m 

To optimize the objective function in (5) with respect to the 
beamwidth (p, let us first assume that p is a constant, and try to 
derive the optimum value of r m . For brevity, in the following 
discussion, we write the objective function as E. Setting the 
derivative with respect to r m as 0, after some calculations we 
have 

dT(\,kr 2 m ) 



Y\-,kriW-p)ipr m + 



= 0, (6) 



where W|,&r^) is calculated as 

r(|, kri) = r [l) + ex P {- kf i) ~ ^ er f ( ■ (7) 



- upper bound 

- simulation 
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Fig. 2. Numerical results and the analytical upper bound for the optimum 
reference distance 



Therefore, 

_^L^^exp(-^)(l-2^ m )- Vfeexp(-<) 

= -2k 3/2 r 2 m exp(-kri). (8) 
Applying (8) to (6), we obtain 

r (! k^j A (1 - p) Wm - 2k 3/2 ri exp (-hi) = 0. (9) 

Since it is difficult to analytically derive the exact solution 
for r m from (9), here we present an upper bound for r m . Since 

T^-,k^>\[T{\,kr 2 m ) + T(2,kr 2 m )} 

= ^p(-kr 2 m )(2 + kr 2 m ), 
using (10) in (9), we have 

A(l - p)k<pr 2 m - 4k 3/2 r m + A(l - p)<p > 0. 

Therefore, 

2k 3/2 - ^4k 3 - 2k[A(l - p)cp] 2 



r m ^ 



kA(l - p)tp 



(10) 



(11) 



(12) 



In Fig. 2, we compare the upper bound of the optimum r m 
with the numerical results for different beamwidth (p, when 
p = 0A. 

B. Jointly Optimizing the Reference Distance and the Trans- 
mission Probability 

Now, let us maximize the objective function by jointly 
optimizing r m and p, for different beamwidth (p. 
Rewrite (5) as 

(3 



E = A 2 p(\-p) exp (~*ptrl) T L,h%\ exp (kr 2 m ) k~ 3/2 sin (|) . 
For brevity, we denote exp(^) as e and Y^kr^A as Y. 



With partial derivatives, we have 
dE 



_ = ^(l-,)sin(|)r^exp(-^) 

w de dY 

• -*x ptrm Ye + T—+e— = 0, (13) 
7i dr m dr m 

This holds only if 



de dY (fi 
Y— + e— = ^Aptr m Ye. (14) 
dr m dr^ n 



Since k = ^^-p + there is f p = To simplify 

things, we can then calculate the derivative with respect to k 
instead of p as 

dE 
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Fig. 3. The optimum transmission probability p vs. the beamwidth 



tr 

1 1 tr, 



2n 



V, l ~ 2 P 

t-n 2 p(l - p) <pA(t - n) 



0. 



(15) 



By using the relationship \r m j£- = fcff , jr m §- = and 
(14), we have 

de dY r m l de dY \ (p 2 
Y— + e— - — \Y- — + e- — = ——Aptr; n Ye. (16) 
dk dk 2k \ dr m dr m J 2nk 



Applying (16) to (15), the following holds: 
trl 3,_i l-2p 2tt 



- V> + l ~ 2p 

t-n 2 



p(l -p) (pMt-n) 

After some calculation, we have 
(ft 



Ye+£-A P triYe = 0. 

(17) 



2 _3 [ J p + 7r/q-7r)](l -2^) 
2(f-*) " m 2 Ml-P) 



(18) 



Given beamwidth ^, we use tp and the optimum to express 
the optimum r m as 



r m = 



1 2[^-7r) + 7T](l-2/7) 3(t-7T) 



^fiX 



p(l-p)t 



(19) 



Interestingly, by jointly optimizing r m and p, we find that 
the optimum transmission probability p in (19) is a constant 
independent of beamwidth (p, which is different from the result 
we had in [1]. The proof is shown in Appendix. 

Since in (19) p is an constant, r m is only related to the 
beamwidth (p. Thus for a given beamwidth (p, there is an 
associated optimum reference distance r m for the transmitters 
to select next hop relay node, as shown in Fig. 4. Also note 
that in (19) r m scales as /i~ 1/2 , which intuitively makes sense. 
This is because as the node density increases, the interferers' 
relative distance to the receiver decreases as VI, it requires a 
shorter transmission distance by the same amount to keep the 
required SINR. By applying (19) in (5), we observe that (5) 
becomes N VI, where N is a constant independent of A. This 
means that the maximum expected density of progress scales 
as @( VI), which conforms to the results in Q and l6l . 
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Fig. 4. The optimum reference distance vs. the beamwidth <p 



IV. Numerical Results and Interpretations 

In this section, we present some numerical results based on 
the analysis in Section III. We choose the path-loss exponent 
a as 3, the node density A as 1, and the outage threshold J3 as 
10 dB. 

In Fig. 3, we show the optimum transmission probabil- 
ity obtained numerically vs. the beamwidth (p. As shown 
in the figure, the optimum transmission probability p is a 
constant, which does not change with the beamwidth (p, this 
confirms our proof in Appendix. Thus, with the selection 
region based routing with directional antennas, no matter 
how much the directional antenna's beamwidth we select, the 
optimum transmission probability always keeps the same as 
p = 0.12. However, when we use the routing strategy with 
omnidirectional antenna, the optimum transmission probability 
p changes with the selection angle (p (see Fig. 5 in [1]). 

In Fig. 4, we compare the optimum reference distance r m 
obtained numerically with that derived in (19), where p is 
chosen optimally as the constant shown in Fig. 3. We see 
that increment of the beamwidth (p leads to the decease of the 
optimum reference distance. This can be explained as follows: 
Increment of the beamwidth (p means more interference seen 
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Fig. 5. Comparison of the expected density of progress for the routing 
protocol using directional antennas with that using omnidirectional antennas 



by a receiver; therefore, the optimum reference transmission 
distance should be decreased to guarantee the quality of the 
received signal and the transmission successful probability. 

In Fig. \5\ we compare the expected density of progress 
of the routing protocol using directional antennas with that 
using omnidirectional antennas in __. From the figure, we 
see that the expected density of progress for selection region 
based routing with directional antennas owns a great advantage 
to that of the routing strategy with omnidirectional antennas. 
This is because directional antennas can bring benefits such 
as reduced interference and increased spatial reuse compared 
with omnidirectional antennas. 

V. Conclusions 

We propose a selection region based multihop routing pro- 
tocol for wireless ad hoc networks with directional antennas, 
where the selection region is defined by the beamwidth and 
the reference distance. By maximizing the expected density 
of progress, we present some analytical results on how to 
refine the selection region and the transmission probability. 
Compared with the routing strategy with omnidirectional an- 
tennas in 0, the routing protocol in this paper owns two 
main advantages: higher expected density of progress and less 
computational complexity when choosing the next hop relay. 

Appendix 

Here, we prove that when jointly optimizing the reference 
distance r m and the transmission probability p, the opti- 
mum transmission probability p is a constant independent of 
beamwidth (p. 

Step 1: We consider the partial derivitive to the reference 
distance r m , which is shown in (9), here rewrite as 

T (| kA A{\-p)<pr m - 2k 3l2 r 2 m exp = 0, (20) 
where k = ^ (f + (1 - p)). Applying k to (20), we have 

fc£exp(-fcy = 0. (21) 



Step 2: Now let us focus on the partial derivitive to the 
transmission probability p. As mentioned in Section III-B, we 
calculate the derivitive with respect to k instead, here rewrite 
(15) as 

dE o (p 
-=^sin(f) 



tr 2 3 

t-7l 2 



1-2/7 



In 



p(l-p) <pA(t-7T) 



(22) 



By using the same notation, defined in Section III-B, i.e., 
e - exp (krfy and T = r (|, krfy, we have 



de 22 
— = exp(£r m )r m . 

Applying (23) and (24) to (22), we have 

tkri 3 1 - 2p pt + 7t(1 - p) 



(23) 
(24) 



3 

'T^~ 2 + p(l-p) 



Ye 



(t-7C) 

+ rexp(^)^-(^) ! =0. (25) 

Step 3: Applying (19) to (21) and (25), we see that in kr 2 m , 
(f can be cancelled, thus (21) and (25) become two equations 
that are independent of (p. Therefore, when jointly optimizing 
r m and p, the optimum p is a constant independent of the 
beamwidth (p, only related to t which is defined in (3a). And 
the numerical value for the optimum p is given in Fig. 3. 
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